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Abstract 

The vacuum-polarization contributions of to the relation between the pole-mass 

M of a quark and the MS parameter rh(M) are evaluated by a straightforward method. 
They are found to approximate very well the exact answer, known through (D(a 2 s ) } thus 
providing a simple physical interpretation. Results are also given for the cases when the 
vacuum-polarization contributions are defined by the pinch technique prescription and spe- 
cific background field gauges. Assuming that the terms n > 3 are also dominant, we eval- 
uate M t /m t (M t ) } compare the results with those of optimization methods, briefly discuss 
Mb/rhb(Mb) } and estimate the irreducible errors in the perturbative series. Implications for 
the electroweak amplitude Ap are emphasized. An update of the QCD corrections to this 
amplitude, including an estimate of the theoretical error, is given. 



1 Leading Vacuum- Polarization Contributions to M/m(M) 

The relation between the pole and running quark masses is a question of considerable practical 
and conceptual interest. For example, in electroweak physics it has been recently observed 
[1-2] that the QCD corrections to the dominant m t -dependent contributions are very small 
through 0(a 2 s ) when the amplitudes are expressed in terms of m t (M t ) (henceforth M and 



rh denote the pole and the MS running masses, respectively). On the other hand, the 
pole mass is the concept most closely related to mass measurements based on kinematical 
considerations [3]. The pole mass is also employed in many calculations of the semileptonic 
and Higss-boson decays, and in the evaluation of rh(M) for the b and c quarks [4]. Recently, 
there has also been considerable interest in studies of contributions to the pole mass associated 
with renormalons, leading to the conclusion that there is an irreducible uncertainty ~ Aqcd 
in the perturbative definition of this concept [5-11]. 

For the ratio M/m(M) there is an important exact perturbative expansion through 0(a 2 s ) } 
due to Gray, Broadhurst, Grafe, and Schilcher [12] : 

_i£_ = 1 + i a(M1+ AV(Af), (1) 

where a(/i) = a s (/i)/7r and K is given by [12,13] 

4 nf fM \ 

K = 16.0065 - n } 1.0414 + 0.1036 + ~J2 A (t7) • ( 2 ) 

3 i=i \M J 

The first term in Eq.(2) corresponds to the quenched approximation, while the second and 
the third are the vacuum polarization contributions of nj massless quarks and the quark of 
mass M, respectively. The A terms represent the mass corrections associated with nj "light 
quarks" with masses M; < M (A(r)/r is roughly constant, being equal to ^ at r = 0, ~ 1.04 
at r 0.3 and n ~ 3 at r = 1). Numerically, K is quite large : 10.96 for the top, 12.5 for the 
bottom, and 13.2 for the charm quarks, respectively. 

Recently [1-2], Eq.(l) in the top quark case has been studied using the BLM [14], PMS 



[15], and FAC [16] optimization methods, neglecting the small A terms. In particular, for 
the BLM procedure, one has 



with jj* t = 0.0960M t . When applied to Eq.(l) in the top quark case, the three optimization 
methods give remarkably close answers, but at the same time they strongly suggest that the 
higher order terms in Eq.(l) have large and increasing coefficients, ~ 10 2 and ~ 10 3 in the 
0(a 3 ) and 0(a 4 ) contributions, respectively. The most rigorous way to settle this issue would 
be to evaluate exactly the next term in Eq.(l). However, because of the difficulty of exact 
multi-loop calculations, prospects for this appear at present to be remote 1 . The question 
naturally arises as to whether there is a simple way to understand large and potentially 
dominant contributions to Eq.(l) occurring in 0(a 2 ) and higher. Following the QED example, 
a natural idea is to incorporate the large effects of one-loop vacuum-polarization bubbles 
(see Fig.l), which may be interpreted as a renormalization of (fi^) 2 , the bare strong coupling 
constant present in the basic one-loop contribution. However, unlike the QED case, the 
concept of vacuum-polarization is ill-defined in non-abelian theories. The reason is that the 
familiar self-energy loops, associated with vacuum polarization, are gauge dependent. To 
underscore the magnitude of this problem, we recall that in the R(_ gauges the logarithmic 
contributions to the one-loop diagrams do not coincide with those associated with the running 
of « s . A significant improvement can be achieved by identifying the one-loop diagrams 
in Fig.l with the pinch technique (PT) self-energy H(k 2 ) [17]. The PT is an algorithm 
that automatically re-arranges one-loop R(_ amplitudes into ^-independent contributions 
endowed with desirable theoretical features. In particular, the logarithms associated with 
the running of a S} as well as the full fermionic contribution, are automatically included in 
H(k 2 ). Defining the self-energy tensor as —i times the associated Feynman diagrams, writing 




(3) 



m t (M t ) 



1 Private communication from D.J. Broadhurst. 
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^-nu(k) = (^g^ — k fl k l/ )H(k 2 ) , and evaluating the PT self-energy H(k 2 ) in dimensional 



regularization, we have 

fl(A; 2 ) = 6( 5s °) 2 (-A; 2 )- e (4 7 r) e - 2 r( e )B(2 - e,2 - e) 



11 -7e 2 ' 



(4) 



1 - e 3 

where k is the external gluon momentum, e = 2 — D/2, -D is the dimension of space-time, and 
B is Euler's Beta function. The first term in Eq.(4) is the ^-independent gluonic contribution 
obtained by adding the pinch parts to the usual self-energy, while the second represents the 
fermionic terms (nj is the number of light quarks regarded as massless). Alternatively, Eq(4) 
can be written in the form 



mk') = 6( S y(-k')-H^f-T(e)BC2 - f,2 - e)b 



1 + 1 ' 



b 1 - ej 

where 6 = 11 — 2nj/3 is the first coefficient of the SU(3) c f3— function. In the e 



using Eq.(8) and implementing the MS renormalization, this becomes 



m 2 )\ 



MS 



«s(M) ( 
47T 



. . ij 2 , 5 4 

ln( -F l + I + I 



(5a) 
limit, 

(5b) 



The mechanism whereby pinching on the relevant Feynman diagrams generates the lowest- 
order H(k 2 ) contribution to the quark self-energy has been shown in Ref. [18]. It has been 
recently pointed out [19] that the PT self-energy coincides with the background-field-gauge 
(BFG) self-energy n(& 2 , £q)bfg evaluated for = 1, where is the gauge parameter 
associated with the quantum loops. At the one-loop level, and for arbitrary £q, the latter is 
given by 2 



n(£ 2 , £ Q ) BFG = 6(^) 2 (-P)- e (4vr) e - 2 r( e )B(2 - e,2 - e) 

xb{l + (4e/6(l - e)) [1 - (1 - £ g )(7 + £g)(3 - 2e)/16]} 



(6a) 



In the e — > limit, employing again Eq.(8) and the MS renormalization, Eq.(6a) reduces to [20] 



2 An expression equivalent to Eq.(6a) has been also derived by G. Weiglein (private communication). 
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A ..ms a s (it) , f, , u 2 . 5 4 L 3 



1-^(1-^X7 + 6?) 



(6b) 



4vr [ v & 2/ 3 6 
Eqs.(5) and (6) indeed coincide for £g = I, as well as for £g = —7 [20]. From Eq.(6b) 
we see that the BFG self-energy also automatically includes the logarithms associated with 
the running of a s , but the non-logarithmic terms depend now on the gauge parameter £q. 
Thus, the concept of vacuum-polarization is also undefined in the BFG approach. A unique 
one-loop answer emerges, however, if one applies the PT within the BFG amplitudes [20], in 
which case one obtains once more 1~I(A; 2 ). 




Figure 1 

Vacuum polarization contributions to the quark self-energy , involving a chain of one-loop 
diagrams. Each bubble can be identified with the PT or BFG one-loop gluon self-energy (see 
text). 



In the first part of our analysis we retain only contributions of 0(6 n_1 a s n (/i)) (n = 1,2,...) 
to M/m(/i), where the expansion is mathematically defined so that there is no residual 
n j— dependence in the corresponding cofactors. The ab-initio rationale is that b is large 
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and that these contributions are uniquely associated with the chains of one-loop bubbles 
depicted in Fig.l, while terms of 0(b l a s n ) (I < n — 2) also arise from other diagrams. A more 
powerful a-posteriori argument is that, as we will explicitly show, the retained contributions 
approximate very well, through (D(a s 2 ) } the exact result of Eq.(l). Consistently with this 
approximation, we first neglect the (1/6) terms in the last factors of Eqs.(5a) and (6a), in 
which limit these expressions coincide. This approximation is mathematically equivalent to an 
approach sometimes referred to as "naive non-abelianization" [8]. The effect of retaining the 
specific (1/6) contributions in Eq.(5a), and in Eq.(6a) for special £q values, will be discussed 
later. 

The sum of contributions involving n = l,2,....,oo bubbles is a geometric series that 
renormalizes {g° s ) 2 according to 

[9s) i-n(P) [) 

In order to express this amplitude in terms of renormalized parameters, we write 

so that the r.h.s. of Eq.(7) becomes 

(47ra s (/i)/Z 3 )(/i 2 eV47r) e 



[l - 3 (6a s (/i)/27rZ 3 ) {-n 2 e</k 2 f T(e)B(2 - e, 2 - e)] 



(9) 



In the MS renormalization, Z 3 is adjusted to cancell the 1/e singularity in the second term 
of the denominator. Thus, for the diagrams under consideration, 

Z 3 =l + ^. (10) 

Setting —k 2 = k 2 , we note parenthetically that in the e — > limit Eq.(7) reduces to a s (e~ 5 / 6 /«), 
the V-scheme running coupling [14]. 

We now turn our attention to the mass renormalization 

M - m = = M) , (11) 



where m is the bare mass and E the quark self-energy. At the one-loop level we have 

d n k ( ?s °) 2 (2M-^(n-2)) 



.4 

= M) = -i- j ^ n k2 ^ 2 + 2p _ k) 



(12) 



t> =M 



where, following the derivation of Eqs.(l,2) [12], we have neglected width effects. In order 
to absorb the vacuum polarization contributions, the next step is to replace (y°) 2 in Eq.(12) 
by the expression of Eq.(9), expand the geometric series and carry out the n-dimensional 
integration. The result is 



4M 00 



/ ba(ji) \ n f(e,ne) 



\ \ 4e ^3 J 



n 



(13) 



where 



f(e,y) 



(&Y [6e><T(l + e)B(2 - e,2 - e)] 



x 2e^T(l + y)r(l - 2y)(l - y)(l - f )/r(3 -e-y). 



(14) 



Recalling Eq.(lO), we now expand (1/Z 3 ) n in powers of x(/-0 = ba(^)/Ae and re-arrange the 
double summation to read 



AM 00 n ( — l) n ~ l 

Z($ = M) = —J2 (X^yf (n - 1)! £ j-^ f(e, le) . 

n=l 1=1 \ / 



(15) 



Next we substitute M — > m in the overall cofactor, as the difference involves terms of 
0(b l a™) (I < n — 2) which are beyond our approximation. Inserting the result in Eq.(ll), 



performing the MS -renormalization 



where 



m = m(n)Z( m ) 



Z( m ) = 1 + a(/i)— + a 2 (/i) (—f - + — ) + ... 



(16) 



(17) 



and dividing by m(/i) we find 



M 



m 



J (m) 



1 + 



S(^ = M) ' 
M 



(18) 



The counterterms present in Z( m ) are important to cancel e singularities but it is easy to 
see that they only contribute finite terms of 0(b l a™) (I < n — 2). For example, Z\ = 0(1) 
and contributes finite terms of 0(a 2 ) but not 0(ba 2 ); Z 22 and Z 2 \ contain terms of 0(b) and 
contribute finite terms of 0(ba 3 ) but not 0(b 2 a 3 ). Thus, in the finite terms of 0(6 n_1 a n ) the 
e n contained in x n (f i ) are cancelled by e factors from /(e,/e). The / summation in Eq.(15) 
can be carried out by expanding 

n 

f(e,le) = J2c P (e)(ley + 0(e n+1 ) , (19) 

p=0 

interchanging the / and p summations so that 

n / -i\n — l n n / -i\n — l 

Sfet7j! /(e -' e) = S^^S&j! ' P + ° (e " +1) ' (20) 

and using the identity 

71 ( — -\) n - 1 ( — -\) n + 1 



Thus 

71 ( — -\\ n - 1 ( — -\\ n + 1 

£ ]J7 — r7\|/( e ' = M~ ^ + c n (e)e n + . (22) 

;=i L \ n V- n - 

The terms of (D(e n+1 ) give vanishing contributions to Eq.(15) as e — > 0. Eq.(21) can be 
proved by applying the differential operator (x-^j to both sides of the binomial expansion 

and then setting x = 1. Inserting Eq.(22) into Eqs.(15,18), we have 



(i-*r = e 0(-!) y > ( 23 ) 



c n (e) H : 



(24) 



From Eqs. (14,19) it follows that 



(\ f( m ^-2^ 7 1-2^ r(l-2e) 

Co(e) = /(e ' 0)= (^^J (w^J r(i + W - e) ■ (25) 
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Expanding 

f(e,0) = J2fre r , (26) 

r=0 

and neglecting terms of 0(b l a s n ) (I < n — 2), only f n and c n (0) contribute to Eq.(24) as 



0. Thus 



(-l) n+1 

Cn(0) + /„ 



(27) 



where the e singularities have explicitly cancelled. The c n (0) can be obtained by using again 
Eqs. (14,19) : 



f(0,y) 



4 e5/3 yr(i + y)r(i- 2y ) 



(28) 



V M2 y (l-y/2)r(l-y) n=0 
We note that the c n (0) depend on /i 2 /M 2 . It is interesting to observe that, using Eq.(28), 

the Borel transform of the first series in Eq.(29), involving c n (0), can be expressed in closed 

analytic form. In fact, it is given by 



If;(|l)" c „ +l(0) = IWM)^l = I 



i^ e 5/3 



2(1 - u)T{u)T{\ - 2u) 1 

u 



M 2 J T(3-u) 
where u = bt/Air and t is the Borel parameter. This expression exhibits the characteristic 
infrared renormalon pole at u = 1/2, and has been already given by Beneke and Braun in the 
leading 1/nj approximation (cf. Eq.(3.13) of Ref. [7]). We have not found, however, a simple 
closed form for the Borel transform of the second series, involving f n . For our purposes, 
which is to evaluate the coefficients of a n } Eq.(27) suffices. 

In summary, the terms of (9(6 n_1 a s n ) are given by Eq.(27). The coefficients f n and c n (0) 
are obtained form Eqs. (25, 26, 28). Setting /i = M, Eq.(27) can be expressed as 

M 



m 



(M) 



where 



A r , 



1 + Mmg^p.j n ^ 



(-If 

2^ Cn+1 ^ + j^TTy. fn+ \ 



(29) 



(30) 
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The coefficients A n are given in Table 1 up to n = 10. We note that the A n are remarkably 
constant for n > 1 and converge to a value 2.301 very close to e 5 / 6 , a factor associated 
with the residue of the infrared renormalon pole. In particular, A\ = 7r 2 /8 + 71/64. The 
second order term is 2bA\a 2 /3 = 1.562056a 2 . Recalling that 6 contains — 2nj/3, we see that 
the coefficient of nj in the a 2 term is -1.0414, in agreement with Eq.(2). This is a welcome 
check, as our calculation should reproduce this factor exactly. For the top, bottom and 
charm cases our approximate evaluation of the a 2 term leads to coefficients 11.98, 13.02, and 
14.06, respectively, to be compared with 10.90, 11.94, and 12.99 from Eq.(2) when the small 
A corrections are neglected. We see that the n = 0, 1 terms in Eq.(29) approximate very 
well the exact results, thus providing a simple physical interpretation of these important but 
technically complex corrections. 

At this stage we return to Eq.(5a), and examine the effect of retaining the specific 1/6 term 
contained in the last factor. It is easy to see that the only changes are an additional factor 
[1 + 4e/6(l - e)]" 1 in the r.h.s. of Eq.(25) and the replacement e 5 / 3 -> e 5 / 3+4 / 6 in Eq.(28). 
This leads to an expansion of the same form as Eq.(29) with A n — > P n} where the first ten 
coefficients are given in Table 2 in the top quark case (6 = 23/3). We see that the P n are 
asymptotically larger than the A n by e 2 / 6 1.30, with smaller differences for low n values. 
In particular, the 0{a 2 (M)) coefficient in M/m(M) is now 1.56205 6 + 9/4. It contains again 
the exact nj term and, in the top-quark case, amounts to 14.23, which is 19% larger that the 
result 11.98 from Eq.(29) and 30% larger than the value 10.90 from the exact calculation of 
Eqs.(l,2). Thus, the P n expansion provides a rough aprroximation to Eqs.(l,2) through 
(D(a s 2 ) } but it is not nearly as precise as Eq.(29). It is also interesting to inquire whether 
the aprroximation of neglecting the 1/6 terms in Eqs.(5a) or (6a) corresponds exactly to a 
specific choice of gauge in the BFG formulation. We note that the gauges £q = —3 ± 4y^2/3 
reduce the expression between curly brackets in Eq.(6a) to 1 + (8/36)(e 2 /(l — e)). As the 
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cofactor of this expression in Eq.(5a) has only a 1/e singularity, in four dimensions these 
particular choices of correspond, in fact, to that approximation. However, the presence of 
the residual e 2 contributions shows that in n— dimensions there exists no BFG that exactly 
reproduces the "naive non-abelianization" prescription. The retention of this e 2 — term does 
not affect Eq.(28), but it introduces an additional factor [1 + 8e 2 /36(l — e)] -1 in Eq.(25). 
This alters the f r parameters and leads to an expansion of the same form as Eq.(29), with 
slightly modified coefficients. In particular the 0{a 2 (M)) coefficient in M/m(M) becomes 
1.56205 b + 1/3. On the other hand, as the c n (0) are not affected, the asymptotic behaviour 
is the same as in the A n expansion. 

As a further application of Eq.(6a), we note that any choice of gauge with £q > 1 will 
increase the coefficient of the e term in the expression between curly brackets, relative to the 
PT result, and therefore lead to asymptotically larger coefficients. In the following we consider 
three interesting cases in which the opposite occurs, namely £q = (which may be referred to 
as the transverse or Landau BFG), £q = —3 and a specific gauge to be defined later. The 
£q = —3 gauge corresponds to the minimum of the e cofactor in Eq.(6a) and, therefore, to 
expansion parameters with the smallest asymptotic behaviour. By a curious coincidence, in 
four dimensions the £q = —3 BFG self-energy also coincides with the one-loop correction to 
the interaction potential between two infinitely heavy quarks (see discussion after Eq.(31)). 
The £q gauges are defined so that the 0{a 2 (M)) contribution in Eq.(29) coincides with the 
exact value K (Eq.(2)) when the small A— terms are neglected. In other words, the complete 
0(a 2 (M)) correctios to M/rh(M) is identified as a self-energy contribution. This can be 
readily implemented by noting that the contribution from Eq.(6a) to the a 2 (M) coefficient in 
Eq.(29) is 1.56205 b + (9 - 23c/3)/4, where c = 3(1 - £ Q )(7 + £q)/16, while Eq.(2) with high 
precision can be rewritten as 1.562056 — 1.07249. Thus, c = 1.73347, which corresponds to 
£q ~ —5.599 or —0.401. We also note that in the general £q gauge, the asymptotic behaviour 
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of the corresponding coefficients in Eq.(29) is e 5 / 6 " 1 " 2 ^ c ^ b . 

For £q = 0, the factor between curly brackets in Eq.(6a) reduces to the expression 
1 - [e/6(l - e)][5/4 - 7e/2]. This introduces a modification e 5 / 3 -> e 5 / 3 " 5 / 46 in Eq.(28) and 
an additional factor [1 — (e/6(l — e))(5/4 — 7e/2)] _1 in Eq.(25). The corresponding L n coef- 
ficients are shown in Table 2 for the top quark case. We see that the asymptotic behaviour 
is reduced by a factor e~ 5 / 86 = 0.922 relative to that in Eq.(29). On the other hand, the 
0{a 2 (M)) coefficient in M/rh(M) becomes 1.562056 — 0.26563, which is very close to the 
result in Eq.(29). 

For general £q, we have the modification e 5 / 3 — > e 5 / 3 + 4 ( 1 - c )/ 6 m Eq.(28) and an additional 
factor {1 + [4e/6(l — e)] [1 — c(l — 2e/3)]} _1 in Eq.(25). As mentioned before, in the case 
c = 1.73347 and we obtain the B* coefficients shown in Table 2 for the top quark. The 
asymptotic behaviour is now reduced by e -1 - 467 / 6 = 0.826 relative to that in Eq.(29). The 
prescription has the interesting feature of incorporating into Eq.(29) the complete 0{a 2 (M)) 
contribution of Eq.(2), which is manifestly gauge invariant. 

For £q = —3, the expression between curly brackets in Eq.(6a) simplifies to 1 — 8e/b. 
Recalling Eq.(8), neglecting the higher-order 1/Z 3 correction and taking the e — > limit, we 
have 



n(^,-3)| BFG = -^6 



. /i 2 5 8 
ln( -p )+ 3"6 



where the superscript MS means that the MS renormalization has been implemented. Eq.(31) 
coincides with the one-loop correction to the interaction potential between two infinitely mas- 
sive quarks [21], and thus it is amenable to direct physical interpretation. The 1 — 8e/b cor- 
rection induces a modification e 5 / 3 — > e 5 / 3-8 / 6 in Eq.(28) and an additional factor (1 — 8e/6) _1 
in Eq.(25). The corresponding V n coefficients are shown in Table 5. Asymptotically, V n ~ 
e~ A l b A n = 0.59A n in the top quark gnificant difference. The 0{a 2 (M)) coefficient in 

M /rh(M) becomes 1.562056 — 7/2, which in the top quark case is 22% smaller than the exact 
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result of Eq.(2). 

In summary, if the basic one-loop vacuum polarization diagram is defined on the basis 
of Eqs.(5a) or (6a) with the neglect of the 1/6 terms, the contributions of (9(6 n_1 a s n ) to 
M/rh(M) are given by the A n expansion of Eq.(29), which corresponds to the so called 
"naive non-abelianization" prescription. If the 1/6 terms in Eq.(5a) are retained, in accor- 
dance with a definition of the basic one-loop vacuum-polarization function based on the 
complete PT self-energy , we have the alternative expansion involving the P n coefficients . 
They are (19-30)% larger than their A n counterparts. If instead, the 1/6 terms in Eq.(6a) 
are retained, corresponding to the BFG framework, the answer depends on the choice of £q. 
We have discussed five particular cases : i) £q = 1 which also leads to the P n expansion, ii) 
£q = —3 ± 4y^2/3, which is very close to the A n case, iii) £q = (Landau or transverse BFG) 
with coefficients L n similar to the A ni but asymptotically reduced by e~ 5 / 86 = 0.922. v) 
which incorporates the exact result of 0{a 2 (M)). iv) £q = —3 which leads to expansion coef- 
ficients V n with the smallest asymptotic limit. In all cases we have the important constraint 
that the a 2 (M) coefficients should be close to the exact answer of Eq.(2). According to this 
criterion, the BFG calculation is singled out, as it reproduces Eq.(2) exactly. However, the 
A n and L n expansions also fare very well, with 0{a 2 (M t )) coefficients that are only 9.8% and 
7.4% larger, respectively, than the exact result. On the other hand, the P n and V n expansions 
are not so accurate, with 0{a 2 (M t )) coefficients 30% higher, and 22% lower, respectively. As 
there are other contributions of 0(b l a s n ) (I < n — 2) which we have neglected, it may be 
argued that the A n expansion represents a more self-consistent approach than its P n} L n and 
V n counterparts. On the other hand, the BFG prescription has the advantage of exactly 
incorporating the most important subleading contributions, namely those of 0{a 2 (M)). It is 
clear that the procedure can be generalized to other QCD calculations where the a 2 (M) 
coefficients are known exactly. 
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After completing our calculation, we have learned that very recently [8] Beneke and Braun 
have also considered the detailed resummation of fermion bubble chains ( characterized as the 
leading corrections in the 1/nj expansion ) and proposed to include the gluonic contribution 
by the heuristic procedure nj — > —36/2, which is refered to as "naive non-abelianization" . In 
the M /rh(M) calculation, this approach should be equivalent to the A n expansion of Eq.(29). 
Their result is expressed in the form 

M _ a(M) ~ fba(M)\ n 

where the large r n coefficients are evaluated by ingenious Borel transform techniques and 
given numerically up to n = 8. We have verified that the relation r n = 2 n+2 n\A n} necessary 
for the consistency of the results of Ref. [8] with Eq.(29), holds numerically with high precision 
up to n = 8. 

Assuming that the terms n > 2 also give dominant contributions, we employ Eq.(29) 
to estimate the higher order coefficients and evaluate M/rh(M). For the top quark case, 
according to Eq.(29) the 0(a 3 ) and 0(a 4 ) coefficients are 86 and 1031, which are close to the 
values 104 and 1041 estimated by optimization methods [1,2]. Their sizeable magnitude 
arises because of the factorial growth and the fact that the expansion parameter involves the 
rather large factor 6/2. As all the terms are of the same sign, Eq.(29) with asymptotically 
constant A n is not Borel-summable. Therefore, the sum is carried out to the optimal term 
of smallest magnitude, which at the same time is regarded as an irreducible theoretical error 
associated with the presumed asymptotic nature of the perturbative QCD expansion [5-11]. 
Although the 0(a 2 ) term in Eq.(29) is a very good approximation, in order to obtain a more 
accurate answer we replace it, in the top-quark case, by 10.90a 2 , the 0(a 2 ) contribution 
from the exact Eq.(2) when the small A corrections are neglected. As an illustration, we 
consider M t = 175 GeV. Using a three-loop f3 function with 5 active flavours normalized to 
a s (mz) = 0.118, we have a s (175GeV) = 0.1074. In this example, the smallest term in Eq.(29) 
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occurs at n = 7 and we obtain M t /m t (M t )= 1.06518 if this term is included, and 1.06482 if 
the series is terminated at n = 6. This is to be compared with the values 1.06470, 1.06470, 
and 1.06463 derived using the three optimization methods [1,2], and 1.05835 from Eq.(l). 
Thus, the QCD correction in the present calculation (0.06482-0.06518) is very close to the 
results found by optimization methods, and differs from the results from Eq.(l) by (11-12)%. 
Table 3 extends this comparison to the range 130 GeV < m t < 220GeV. The second column 
gives the value from Eq.(l), the third column that from the BLM expression (Eq.(3)), the 
fourth and fifth columns give the results of incorporating the higher order terms in Eq.(29) 
up to n = 6 and n = 7 (previous to smallest and smallest terms, respectively). The smallest 
term for M t = 175GeV is 3.52 xl0~ 4 , so that our estimate of the irreducible theoretical error 
in M t is 3.52 X 10~ 4 X 175GeV/1.065 ~ 58MeV. If we use Stirling's approximation and the 
one-loop relation 2/ba(M) = ln(M/AQCD), our irreducible error estimate can be written as 

6MK 5(-f)) 1/2 e»/«A QCD (32, 

In this case Aqcd~ 85 MeV (from the one-loop relation) and from (32) one has SM t 62 
MeV. An alternative procedure to estimate SM is based on the ambiguity of the Borel trans- 
form. For example, using Eq.(34) of Ref. [8], one obtains an expression similar to Eq.(32) 

l/2 

with (a s (M)/b) — > 1/6. In the top quark case this increases the uncertainty by a factor 
(a s (M) b) 1.10 so that SM 68 MeV. Fortunately, such uncertainties are phenomeno- 

logically negligible. Because complete higher order results are not available, the actual error 
in current calculations is, of course, likely to be significantly higher. In Ref. [2] the last term 
in Eq.(3) was used as an estimate of the magnitude of the theoretical uncertainty and this 
gives ±2.74 X 10~ 3 . It is worth noting that an analogous use of the P n and V n expansions, 
i.e. the two extreme cases we have considered, lead to values of M t /m t (M t ) that differ from 
the BLM result (Eq.(3)) by (2 — 2.5) X 10~ 3 . On the other hand, the L n and B* n expansions 
differ from the BLM result by ~ 1 X 10~ 3 . We also emphasize that the optimization methods 
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applied directly to the 0(a 2 s ) expansion (Eq.(l)) fail to uncover the factorial growth displayed 
in Eq.(29). For instance, if in Eq.(3) is expanded in terms of a(M t ), it obviously gen- 

erates a geometric series. The numerical closeness of the two calculations is therefore very 
interesting. 

We briefly illustrate the application of Eq.(l) and Eq.(29) to the bottom quark case using 
Mf, = 4.72 GeV [13], in which case a s (4.72GeV) = 0.2162. We again improve the accuracy by 
replacing the a 2 term in Eq.(29) by 12.5a 2 , the result obtained from the exact Eq.(2) including 
the A corrections. In this case, the smallest term is n = 3, and we have Mb/rhb(Mb) = 1.214 
if it is included, and 1.184 if the sum is stopped at n = 2. This is to be compared with 1.151 
from Eq.(l). Such differences are not too surprising. We note that the expansion from Eq.(l) 
is 1 + 0.0918 + 0.0592. Thus, although Eq.(l) is routinely applied to the b and c quarks [4], 
its convergence properties suggests a large theoretical uncertainty in those cases. As the value 
of rhb(Mb) is frequently derived from M& by using Eq.(l) [4], one obtains in this case rhb{Mb) 
= 4.72 GeV/1.151 = 4.10 GeV, while application of the present n = 2 calculation leads to 
m b (M b ) = 4.72 GeV/1.184 = 3.99 GeV, a difference of 114 MeV or « 3%. The smallest 
term in the series amounts in this case to 2.97xl0~ 2 , so that our estimate of the irreducible 
theoretical uncertainty in M b is 2.97 xl0~ 2 x 3.99 GeV ~ 119 MeV, which is twice as large as 
that given in Ref. [6]. This difference can be roughly traced to the factor e 5 / 6 in Eq.(32). A 
more detailed treatment of M/m(M) for the b and c cases has recently been proposed by 
Beneke and Braun [8]. In the b case their central value is close to the average of the two 
values we have given, while their irreducible error is about 25% smaller. 

2 QCD corrections to (Ap)f 

The study of the QCD corrections to (A/?)/, the fermionic component of A/9, is a subject 
of considerable interest, as this amplitude contains the leading asymptotic contributions, 
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for large M t} of the basic radiative corrections Ar,Ar, and A/3 [22,23]. We also recall that 
pj = [1 — (A/?)/] -1 is frequently separated, as an overall renormalization factor, when neutral 
current amplitudes are expressed in terms of G^. In this section we apply the previous results 
to the study of these corrections. We also update the discussion of Refs.[l,2], including an 
estimate of the theoretical error, to take into account very recent modifications in the (D(aa 2 s ) 
calculations. 

Neglecting higher-order electroweak effects ~ 0(a/Trs 2 ) 2 (M 2 /m 2 v ) } but retaining QCD 
corrections, (A/?)/ can be written as 

(Apl/ = S [1 + * Q0Dl • (331 

where the first factor is the one-loop result and Sqqd represents the relevant QCD correction. 



Alternatively, calling p t the solution of rh(p) = /i, where rh t {p) is the MS -running mass, 
Eq.(33) can be expressed in the equivalent forms 



8V2tt 
and 



( = 8 ^ 2 [ + QCD] ' ( } 



Eq.(34) can be interpreted as a pure-MS expression, while Eq.(35) is very useful because it 
can be used in conjunction with the results of Section 1. 

Two recent evaluations of the complete three-loop corrections of (D(aa 2 s ) have been given 
by Avdeev, Fleischer, Mikhailov, and Tarasov [24], and by Chetyrkin, Kiihn, and Steinhauser 
[25], in the limit M& = 0. After an initial discrepancy the two results now agree and, to good 
accuracy, are given by 

6 QCD = -2.8599 a(M t ) - (23.525 - 1.7862 n s )a 2 (M t ) , (36) 

where the 0{a(M t )) term is the Djouadi-Verzegnassi result, and nj = 5 is the number of 
massless flavours. Numerically, the second order coefficient is —14.594. Using Eq.(l) in the 
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A = approximation and the relation fi t = [1 + (8/3)a 2 (M t ) + ...]m(M t ), Eq.(36) implies 

^ C 5 D = -0.19325 a(M t ) - 3.970 a 2 (M t ) , (37) 

Aqcd = -0.19325 a(M t ) + 1.364 a 2 (M t ) . (38) 

As pointed out in Ref. [25], the second order coefficient in Eq.(37) is almost entirely due to 
the opening of a new channel in (D(aa 2 s ) } namely the double triangle graph. It is interesting 
to compare the new results of Eqs. (36-38) with a very simple estimate made before the 
complete three-loop calculations were carried out [26]. In that analysis the correction Aqcd 
was written as 

Aqcd = -0.19325 a(M t ) + C a 2 (M t ) , (39) 

and the constant C estimated to be C = ± 6 on the basis of convergence assumptions 
(optimization arguments carried out in Ref. [26] also led to C ~ +3, but the more conservative 
estimate ± 6 was employed in the analysis). In turn, C = ± 6 implies a second-order 
coefficient —15.958 ± 6 in Eq.(36). Thus, we see that the central values in the estimate of 
Ref. [26], namely C=0 in Eq.(39) and —15.958 in Eq.(36), are amusingly close to the new 
revised second-order coefficients in the exact calculation . 

Applying directly the optimization methods to Eq.(36), we find 

«5 QCD = -2.8599 a(0.1536M t ) + 5.947 a 2 (0.1536M t ) (BLM) , (40) 

«5 QCD = -2.8599 a(0.2241M t ) + 1.803 a 2 (0.2241M t ) (PMS) , (41) 

«5 QCD = -2.8599 a(0.2642M t ) (FAC) . (42) 

Eqs. (40-42) can be obtained, for example, by using the explicit expressions given in Eqs. (16- 
18) of Ref. [2]. As an illustration, we consider the case M t = 175 GeV. Employing a three-loop 
f3 function with five active flavours ( we neglect the very small three-loop discontinuity at 
fi = Mt) and normalizing once more a s (fi) such that a s (Mz) = 0.118 for Mz = 91.19 GeV, 
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Eqs. (40-42) give, for M t = 175 GeV, <5 QCD = -0.1192, -0.1198, and -0.1197, respectively. 
This is to be compared with —0.1149 from Eq.(36). It is interesting to note that: i) The BLM 
approach works considerably better with the recently modified version of <5qcd, in the sense 
that the residual 0{a 2 ) term is not as large as before, ii) The three optimization methods 
give now close answers, iii) The difference between the results from the optimization methods 
and Eq.(36) is 4.9 X 10~ 3 , which is roughly as large as the effect of the recent modifications, 
namely 4.04 a 2 (M t ) = 4.7 X 10~ 3 . It should be stressed that the 4.9 X 10~ 3 variation is due 
to different ways of accounting for the contributions of 0(a 3 ) and higher. 

The presence of large and increasing coefficients, the observations in iii), and the large 
vacuum-polarization contributions to M t /m(M t ) we have encountered in Section 1, strongly 
suggest that the corrections of 0(a 3 ) and higher to Eq.(36) are significant. It is therefore a 
good idea to examine alternative evaluations of Sqqd that involve second order coefficients 
of C(l), rather than 0(10). Such strategy was explained in Refs.[l,2], and is updated in the 
following to take into account the very recent modifications in the 0{aa 2 s ) corrections and 
the results of Section 1. 

Combining Eqs. (33) and (35) we have 

1 + ^QCD = (^^) [1 + A Qcd] • (43) 

The correction Aqcd is very small. However, in order to incorporate large vacuum-polarization 
effects that are induced in (9(a 3 ), we replace Eq.(38) by the more detailed expression 

16 

Aqcd = -0.19325 a(M t ) - 3.970 a 2 (M t ) + — a 2 (0.252 M t ) - 9.97 a 3 (0.252 M t ) . (44) 

This is obtained by combining Eq.(37) with the BLM optimized expansion 

—7^- = 1 + ^a 2 (0.252 M t ) - 4.47a 3 (0.252 M t ) , (45) 
m t {M t ) 3 

derived in Ref. [1]. The large vacuum-polarization effects of 0(a 3 ) incorporated in Eq.(45) 
are due to the presence of the pole mass in the argument of rh. As Eq.(45) is exactly known 
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through 0(a 3 ) [1], this is a very instructive illustration of the fact that the occurence of the 
pole mass M in theoretical expressions is very likely to engender large vacuum-polarization 
effects, as well as large coefficients in higher orders, if a s (M) is employed as an expansion 
parameter. The correction Aqcd, evaluated on the basis of Eq.(44), is given as a function 
of M t} in Table 4. We see that Aqcd ~ —(2.3 — 2.6) x 10~ 3 is very small. This is close in 
magnitude to the values reported before [1], but of opposite sign. The difference is due to 
the modification of the second-order coefficients in the exact calculation . 

Table 5 compares several evaluations of <5qcd, as a function of M t . The second column 
gives the values from Eq.(36) while the remaining ones report the results from using Eq.(43), 
Table 4 for Aqcd, and the evaluations of M t /m t (M t ) given in Section 1. Specifically, the third, 
fourth, and fifth column employ the M t /m t (M t ) values from the BLM approach (Eq.(3)), and 
the n = 6 and n = 7 summations explained in Section 1, respectively (see Table 3). It is 
clear that the last three columns, and especially the third and fourth, are quite close. On 
the other hand, they are larger than the first by roughly (4-6)% depending on M t and the 
particular entries being compared. This is not too surprising because the last three columns 
incorporate relatively large vacuum-polarization contributions of 0(a 3 ...) to M t /m t (M t ) } not 
included in the first one. Because of the closeness of the last three columns, there is no point 
in drawing distinctions among them and, for definiteness, we adopt the entries in the third 
one, based on Eq.(43) and the BLM optimization of Eq.(3), as our central values. Following 
Ref. [2], we estimate the theoretical error due to unknown higher order corrections by the 
magnitude of the last retained terms. We emphasize that, although the three optimizations 
of M t / m t (M t )give nearly the same numerical answer, in order to analyze the error we employ 
the BLM approach, which leads to the most conservative estimate. From Eq.(3) we have 

6 (^] « ± 2 , Xl ;° 73a2 ^ ) - ±1-78 «V) • (46) 
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From Eq.(38) 

«5(A QCD ) ~ ±10 a 3 (0.252 M t ) . (47) 
Inserting these uncertainties in Eq.(43) and adding them linearly we obtain 

S(S QC d) = ±(4.8, 4.9, 5.2, 5.9) x 10" 3 , (48) 

for M t = 220, 200, 175, 130 GeV, respectively. These error estimates coincide with the 
values reported in Ref. [2]. 

In analogy with the discussion in Refs.[l,2], the values for Sqqd given in Table 5 can 
accurately be represented by simple empirical formulae. We find that our central values (third 
column of Table 5) and error estimates due to higher order corrections can be conveniently 
expressed as 

S QC D = -2.8599 a{iM t ) , (49) 

i = 0.2601°;™ , (50) 

while Eq.(36) (second column of Table 5) corresponds to £ = 0.339. We emphasize that 
Eqs. (47,48) are not the result of a FAC optimization. They are simply empirical formulae 
that reproduce the values in the tables with errors of at most 1 X 10~ 4 for £ = 0.260 and 
£ = 0.339, and at most 2 x 10" 4 for £ = 0.204. 

Combining quadratically the theoretical error in Eq.(48) with that arising from the ±0.006 
uncertainty is a s , we obtain an overall uncertainty in Sqqd of ±(8.0, 8.2, 8.6, 9.7) X 10~ 3 for 
Mt = (220, 200, 175, 130)GeV. The effect on the predicted values of mw and sin 2 0w(mz) 
due to a shift in Sqqd can be obtained using Eqs. (37a) and (37b) of Ref. [27]. We have 

8m w c 2 x t 8(8qcD) 



mw 2(c 2 — s 2 — 2c 2 x t ) 



(51) 



6s 2 = --^^-^^ls(S QCD ) , (52) 
167T c 2 — s 2 m 2 z 
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where x t = 3G |U M t 2 /(8v^7r 2 ), s 2 = 1 — c 2 = 1 — m^/m|, and s 2 = 1 — c 2 = sin 2 9w(mz) and 



a = a(mz) ~ 1/127.9 are the MS parameters. For the above M t entries and Mjj = 300GeV, 
and using the overall uncertainty <5(<5qcd), we obtain Smw = ±(7.1, 6.0, 4.8, 3.0) MeV and 
Ss 2 = ±(4.0, 3.4, 2.8, 1.7) X 10~ 5 , with a very similar shift for the effective electroweak 
parameter sin 2 #^. The effect on the M t prediction from precision electroweak physics can 
be inferred from detailed tables such as those given in Ref. [27]. However, a simple estimate 
can be obtained by noting that, to a good approximation, the precision electroweak data fixes 
the value of (Ap)j. Recalling Eq.(33), we have therefore the approximate relation 

gMt ^ .^J^L , (53) 

2 [1 + 6qcd\ 

from which we find 6M t « ±(9.9, 9.3, 8.6, 7.2) x 10" 1 GeV. 

It is interesting to note that: i) The results of the a(M t ) expansion (Eq.(36)) lie almost ex- 
actly at the lower boundary (in absolute value) of the band derived by our central values (third 
column of Table 5) and error estimates (Eqs. (48-50)). ii) On the other hand, the values from 
the optimized versions of Eq.(36), namely Eqs. (40-42), are not only close among themselves, 
but also lie near the central values reported in this paper, iii) The latter represent a fractional 
enhancement of (24-22)%, depending on M t} over the one-loop Djouadi-Verzegnassi result, 
iv) Our value for M t = 200 Gev, namely Sqqd = —0.1174 ± 0.0049, is quite close to the 
estimate Sqqd = —0.121 ± 0.006 given in Ref.[25] before the complete 0{aa 2 s ) result became 
available. 
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Table 1 

The A n coefficients in Eq.(29). This expansion corresponds to the neglect of all contributions 
of 0(b l a s n ) (l<n- 2). See text. 



n 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


A n 


1 


2.34308 


2.20283 


2.28873 


2.27584 


2.30200 


2.29632 


2.30185 


2.30003 


2.30127 


2.30077 



Table 2 

Alternative prescriptions for the expansion coefficients in Eq.(29). i) The P n replace the 
A n when the 1/6 terms in the PT self-energy are incorporated, ii) The L n coefficients are 
similarly based on the transverse or Landau (£q = 0) BFG. iii) The B* correspond to the 
£q prescription, iv) The V n correspond to the £q = —3 BFG. Among the BFG coefficients, 
the V n have the smallest asymptotic limit. These alternative prescriptions are nj— dependent 
and are given here for the top quark case (b = 23/3). 



n 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Pn 


1 


2.78329 


2.84487 


2.93323 


2.95309 


2.98043 


2.98150 


2.98642 


2.98592 


2.98687 


2.98664 


L n 


1 


2.29111 


2.02192 


2.11758 


2.09633 


2.12387 


2.11601 


2.12210 


2.11980 


2.12122 


2.12060 


K 


1 


2.13324 


1.79728 


1.91085 


1.87524 


1.90510 


1.89496 


1.90208 


1.89907 


1.90076 


1.89998 




1 


1.65829 


1.23658 


1.41819 


1.32999 


1.38095 


1.35691 


1.36964 


1.36347 


1.36663 


1.36507 
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Table 3 

Comparison of different evaluations of M t /m t (M t ). The second column gives the results of 
Eqs.(l,2), while the third one is based on the BLM optimization Eq.(3). The fourth and fifth 
columns are based on Eq.(29) with the summation truncated at n = 6 and n = 7, respectively, 
and the second order coefficient replaced by its exact value from Eqs.(l,2). 



M t /m t (M t ) 



M t 


Eqs.(l,2) 


Eq.(3) 


Eq. (29) 


Eq. (29) 


(GeV) 






n = 6 


n = 7 


130 


1.06140 


1.06875 


1.06902 


1.06951 


150 


1.05989 


1.06673 


1.06693 


1.06735 


175 


1.05835 


1.06470 


1.06482 


1.06518 


200 


1.05708 


1.06303 


1.06311 


1.06342 


220 


1.05621 


1.06190 


1.06195 


1.06222 
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Table 4 



The correction Aqcd (Eq.(44)) 



M t 
{GeV) 


10 3 A QCD 


130 


-2.27 


150 


-2.39 


175 


-2.50 


200 


-2.59 


220 


-2.64 



Table 5 

Comparison of several evaluations of Sqqd ■ 



<*>QCD 



M t 


Eq. (36) 


Eqs.(3,A3) 


Eqs.(29,A3) 


Eqs.(29,A3) 


{GeV) 






n = 6 


n = 7 


130 


-0.1205 


-0.1265 


-0.1269 


-0.1277 


150 


-0.1177 


-0.1233 


-0.1236 


-0.1243 


175 


-0.1149 


-0.1200 


-0.1202 


-0.1208 


200 


-0.1125 


-0.1174 


-0.1175 


-0.1180 


220 


-0.1109 


-0.1155 


-0.1156 


-0.1161 
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